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We study the correlation length of the two-dimensional Edwards- Anderson Ising spin glass with 
bimodal interactions using a combination of parallel tempering Monte Carlo updates and a rejection- 
free cluster algorithm in order to speed up equilibration. Our results show that the correlation length 
grows ~ exp(2J/r) suggesting through hyperscaling that the degenerate ground state is separated 
from the first excited state by an energy gap ~ 4J, as would naively be expected. 
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I. INTRODUCTION 

Despite intense research, ^-^ spin glasses still pose many 
unanswered questions. Thus it is of great interest to un- 
derstand the nature of the spin-glass state for realistic 
short-range spin glasses. A model that has proven to 
be a good "workhorse" due to its simplicity and ease 
of implementation is the Edwards- Anderson Ising spin 
glasE'^ with bimodal interactions. Here we consider this 
model in two dimensions, for which the spin-glass tran- 
sition occurs at zero temperature. Nonetheless the low- 
temperature properties of the model remain controver- 
sial: There are different predictions on the energy gap 
between the ground state and the first excited state, and 
the critical behavior of the correlation length is not fully 
understood. 

Wang and Swendseni first suggested an exponential 
scaling of the specific heat of the two-dimensional bi- 
modal Ising spin glass: 

Cv - T-^e-^f^^ , (1) 

where (3 = 1/T represents the inverse temperature, J 
is the mean interaction strength of the (random) bonds, 
and ^ is a numerical prefactor. Although they argue an- 
alytically that the energy gap should be 4 J, i.e., A — i, 
their final estimate from their numerical work is ^ = 2. 
In addition, they estimated the correlation length expo- 
nent v to be finite and positive, implying power-law scal- 
ing of the correlation length. Unfortunately, their simu- 
lations were done on small system sizes with few disorder 
realizations, suggesting that their results are governed by 
corrections to scaling. 

Saul and Kardar^.^ by means of an exact integer al- 
gorithm to estimate the partition function of the two- 
dimensional ± J spin glass, later argued that A = A and, 
in addition, proposed an exponential scaling form for the 
correlation length, i.e., ^ ~ exp(n/3J) with n = 2. The 
exponential scaling present in the bimodal glass was later 
verified by Houdayer— by analyzing the scaling properties 
of the spin-glass susceptibility and the Binder cumulant 
for large system sizes using a novel cluster algorithm. Re- 
cently, Lukic et al& reconsidered the problem by comput- 
ing the exact partition function of the system for system 



TABLE I: Different estimates for n and A. From finite- 
size scaling arguments in two dimensions, we expect that the 
specific heat scales as Cv ~ where ^ is the correlation 
length. Following the exponential scaling for Cv suggested 
by Saul and Kardar (Ref. |^, we expect in general Cv ~ 
exp{—Af3J) and ^ ~ exp(n/3J), with A = 2n. 



Reference n A 
Wang and Swendsen (Ref. ^ 2 

Saul and Kardar (Ref. 5) 2 4 

Houdayer (Ref.^ 2 

Lukic et al. (Ref. ^ 12 

Katzgraber and Lee (this work) 2 



sizes up to 50 x 50 spins. They conclude that A = 2. 
According to hyperscaling, the singular part of the free 
energy scales as ^"^ (in two dimensions) and so, if ^ di- 
verges exponentially, the leading exponential dependence 
of the specific heat will also be given by Therefore, 
one expects A ^ 2n, and in the case of Lukic et al. this 
would mean that ^ ~ exp(n/3J) with n = 1. For a sum- 
mary of all predictions of A and n, the prefactors to f3J 
in the exponential scaling of specific heat and correlation 
length, respectively, see Tabled 

In this work we compute the finite-size correlation 
length^' directly via Monte Carlo simulations for large 
system sizes and show that n changes continuously for 
intermediate system sizes, yet saturates at n « 2 for 
large enough systems. This suggests via hyperscaling 
that A = A, i.e., the energy gap in the two-dimensional 
Ising spin glass is ~ 4 J. 

The paper is structured as follows. In Sec. ^ we intro- 
duce the model and observables and in Sec. lIIII we present 
our results. Conclusions are summarized in Sec. IIVI 

II. MODEL AND OBSERVABLES 

The Hamiltonian of the two-dimensional Ising spin 
glass is given by 

H ^ -'^Ji■iS^Sj, (2) 
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TABLE II; Parameters of the simulations. A'^samp represents 



the number of disorder reahzations computed; N^, 



the to- 



tal number of Monte Carlo sweeps of the 2Nt replicas for a 
single sample. Nt is the number of temperatures in the paral- 
lel tempering method and Tmin represents the lowest temper- 
ature simulated. For L < 16 we have used standard parallel 
tempering Monte Carlo simulations, whereas for L > 24 we 
have used a combination of parallel tempering and cluster 
updates (Ref. 0). 



L 


samp 


sweep 




At 


4 


10 000 


2.0 X 10^ 


0.050 


20 


6 


10 000 


2.0 X 10^ 


0.050 


20 


8 


10 000 


2.0 X lO'^ 


0.050 


20 


12 


10 000 


4.0 X 10^ 


0.050 


20 


16 


10 000 


1.0 X 10® 


0.050 


20 


24 


5 003 


2.0 X 10*^ 


0.050 


20 


32 


5 031 


2.0 X 10® 


0.050 


20 


64 


500 


4.2 X 10® 


0.200 


39 


96 


609 


6.5 X 10® 


0.200 


63 


128 


420 


2.0 X 10® 


0.396 


50 



where the sum ranges over nearest neighbors on a square 
lattice with periodic boundary conditions. Si represent 
Ising spins taking values ±1, and the interactions J^- are 
bimodally distributed, i.e., Jij G ±1. For the Monte 
Carlo simulations we use a combination of single-spin 
flips, parallel tempering updateSfiiii^ and rejection-free 
cluster movesSii^ in order to speed up equilibration. To 
ensure that the system is equilibrated, we perform a loga- 
rithmic data binning of all observables (energy, spin-glass 
susceptibility, and correlation length) and require that 
the last three bins agree within error bars and are inde- 
pendent of the number of Monte Carlo sweeps A'swcop- 
The parameters of the simulation are listed in Table Ull 

As mentioned in SecU according to Saul and Kardarji^ 
the correlation length of the two-dimensional ±J spin 
glass scales as 



(3) 



where /3 = 1/T and n = 2. In this paper we determine 
the value of n from Monte Carlo simulations. 

In order to compute the correlation length ^, as well 
as estimate the critical exponent ?7, we compute the spin- 
glass susceptibility 



XSG = N[{q^)] 



(4) 



where [• • -j^^ represents a disorder average and (■ ■ •) a 
thermal average, and 



(5) 



is the Edwards- Anderson spin-glass order parameter. In 
Eq. ^ N = represents the number of spins and a and 



f3 represent two replicas of the system with the same dis- 
order. Due to the exponential scaling^ of the correlation 
length we expect^ 



Xsg{T,L) = L^-^C 



[3- - \nL 
n 



(6) 



The finite-size correlation lengtl&iSiiiiiSiiSiii is given 
by 



2sin(|k,„in|/2) 



Xsg(O) 

XSG(kmin) 



-,1/2 



1 



(7) 



where kmin — (27r/i,0, 0) is the smallest nonzero wave 
vector, and xsG(k) is the wave- vector-dependent spin- 
glass susceptibility: 



(8) 



Because £,l/L is dimensionless, we expect 
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(Tc = 0) . (9) 



and 



Here X is a scaling function. Note that in Eqs. © 
we have explicitly left n as a "variable." 

To determine the asymptotic value of n, we compute 
the bulk correlation length at low temperatures using 
the method of Kim^"' (first introduced in Ref . fisi) . The 
finite-size scaling relation for the correlation length can 
be written as 



Coo p f £,L 



(10) 



where we determine f{x) by fitting to data in the range 
0.500 < T < 1.391 where we have data for the correlation 
length in both the bulk and finite-size regimes. Using 
/(.t), we then determine from Eq. I|10|) using data for 
L = 96 and 128 in the range 0.396 < T < 1.391. 



III. RESULTS 

In Fig. n we show data for the natural logarithm of the 
finite-size correlation length as a function of 1/T for dif- 
ferent system sizes. According to Eq. Q, we expect data 
for In(^L) vs 1/T to asymptotically approach a slope of n. 
The data presented in Fig. ^ show good agreement with 
n « 2, which in turn agrees with the results of Refs. 
and ,6j We have also used the extrapolation method by 
Kim^^ and show data for the bulk correlation length ex- 
trapolated from data for the largest system sizes. These 
data also agree well with the naive scaling with n = 2, 
and not with n = 1 (dotted line in Fig.^, as predicted 
by Lukic et al& from specific heat studies. Corrections 
to scaling in the two-dimensional bimodal spin glass are 
strong for system sizes L up to 64. Hence the results of 
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FIG. 1: (Color online) Data for the finite-size correlation 
length: ln(^_t) vs l/T. The slope of the curves in the thermo- 
dynamic limit determine n. For small values of l/T the data 
are independent of system size L, but the data "peel off" from 
this common curve at a value of 1 /T which increases with in- 
creasing size. The slope of the common curve asymptotically 
seems to approach the value n = 2 (dashed line). This is sup- 
ported by data for the bulk correlation length (•) which shows 
an asymptotic slope also compatible with ^ ~ exp(2/3J), i.e., 
n = 2. (The dotted diagonal line has slope 1.) 



Lukic et al. for L < 50 could be influenced by corrections 
to scaling that mask the true critical behavior. 

In order to test this hypothesis, we show a finite-size 
scaling plot of the data of the finite-size correlation length 
according to Eq. ||5J) with n = 2 (Fig.[21). We see that the 
data for L ^QA are not in the asymptotic limit, whereas 
the data for L > 64 scale reasonably well. This shows 
that in order to obtain precise estimates of n, very large 
system sizes are required]^ In addition, the exponential 
scaling means that, for the two-dimensional bimodal spin 
glass, the critical exponent for the correlation length is 
infinite, i.e., 

v^oo. (11) 

We have also tried to scale the data according to Eq. @ 
using n = 1 in accordance with the prediction of Lukic et 
al.^ As shown in Fig.Othe data do not scale well for any 
system size and any temperature suggesting that n = 1 
is not probable. 

Data for the spin-glass susceptibility show similar 
finite-size effects as in the case of the correlation length. 
In Fig. 0] we show a finite-size scaling plot of the spin- 
glass susceptibility xsG assuming n = 2 and 77 = 0.138. 

Note that scaling works "reasonably well" for a large 
range of 77 values (not shown) thus illustrating again^^ 
that the spin-glass susceptibility is not a good observable 



FIG. 2: (Color online) Parameter-free scaling plot of the 
finite-size correlation length according to Eq. Q using n = 2. 
Only the data for large values of L (L > 64) scale for low 
temperatures. 




FIG. 3: (Color online) Parameter-free scaling plot of the 
finite-size correlation length according to Eq. JUJ using n = 1. 
The data do not scale well for any temperature and system 



with which to study critical properties of glassy systems. 
At T = Tc (= 0) we expect, according to Eq. jnj, 



XSG 



L^"" (T = Tc 



(12) 



In order to estimate the critical exponent 77, we plot 
ln(L~^XSG) vs ln(L) and extract the slope of the 
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IV. CONCLUSIONS 
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FIG. 4: (Color online) Scaling plot of the spin-glass suscepti- 
bility with Tj = 0.138 and n — 2. The inset shows the efltective 
critical exponent rjesiT). 



curves, which should be proportional to a temperature- 
dependent effective exponent ~'qcs{T), where 



hm ■qc'niT) = 77 



(13) 



The inset of Fig. 01 shows the effective exponent rjcs (T) 
as a function of temperature. We estimate from our data 



•q = 0.138 ±0.005 , 



(14) 



a result which is slightly smaller than a recent estimate by 
HoudayeriS We expect the value of the critical exponent 
•q to change slightly with larger system sizes. However, 
because we do not have data for the largest L at the 
lowest temperatures (see Table 1^), we cannot make a 
strong statement regarding 77. 



To conclude, we have shown that the correlation length 
of the two-dimensional Ising spin glass with bimodal in- 
teractions scales exponentially ~ exp(n/3J) with n = 2. 
These results are in agreement with previous work by 
Saul and Kardari and HoudayerA However, if one as- 
sumes the hyperscaling result A = 2n, our results dis- 
agree with specific heat studies by Lukic et al& as well as 
Wang and Swendseni who find A = 2. Using the hyper- 
scaling relation, our results also imply that the excitation 
gap for the bimodal spin glass is ~ AJ . Exponential scal- 
ing of the correlation length means that the critical expo- 
nent for the correlation length is infinite {v — 00) for the 
two-dimensional bimodal spin glass. In addition, we esti- 
mate the critical exponent 77 and find that rj « 0.138(5). 

Although in this work we were able to present data 
that can be extrapolated to the bulk regime for a certain 
temperature range therefore allowing us to draw conclu- 
sions in the thermodynamic limit for the scaling of the 
correlation length, finite-size effects are very strong in 
this system and so an analysis with yet larger system 
sizes at lower temperatures is desirable. While it is un- 
likely, a change in the asymptotic behavior at larger L 
cannot be ruled out completely. 
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